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Modelling the gluon propagator* 
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Scaling of the Landau gauge gluon propagator calculated at /3 — 6.0 and at ft = 6.2 is demonstrated. A variety 
of functional forms for the gluon propagator calculated on a large (32 3 x 64) lattice at /? = 6.0 are investigated. 



1. Scaling behaviour 

In this note we focus on modelling the gluon 
propagator calculated in We use the same 
conventions and definitions here as in Qj. 

The lattice gluon propagator is related to the 
renormalised continuum propagator D R (q; fi) via 
a renormalisation constant, 



D L (qa) = Z 3 (/j,a)D R (q;/j) 



(1) 



The asymptotic behaviour of the renormalised 
gluon propagator in the continuum is given to 
one-loop level by 



D R (q) 



1 /l 



ln(g 2 /A 2 



-dr 



(2) 



r \2 

= 13/44 in Landau gauge for quenched 



with du 
QCD. 

Since the renormalised propagator D R (q;fj,) is 
independent of the lattice spacing, we can use (Q) 
to derive a simple, ^-independent expression for 
the ratio of the unrenormalised lattice gluon prop- 
agators at the same value of q: 



D c {qa c ) Z 3 (n,a c )D R (q;ii)/a 2 c 



Z f a 



(3) 



where the subscript / denotes the finer lattice 
(/3 = 6.2 in this study) and the subscript c de- 
notes the coarser lattice (j3 — 6.0). We can use 
this relation to study directly the scaling proper- 
ties of the lattice gluon propagator by matching 
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the data for the two values of (3. This matching 
can be performed by adjusting the values for the 
ratios Rz — Zf/Z c and R a = a,f/a c until the two 
sets of data lie on the same curve. We have im- 
plemented this by making a linear interpolation 
of the logarithm of the data plotted against the 
logarithm of the momentum for both data sets. 
In this way the scaling of the momentum is ac- 
counted for by shifting the fine lattice data to the 
right by an amount A a as follows 



]nD c (]n(qa c )) = lnD/(ln(ga c ) 



(4) 



Here Az is the amount by which the fine lattice 
data must be shifted up to provide the optimal 
overlap between the two data sets. The matching 
of the two data sets has been performed for values 
of A a separated by a step size of 0.001. Az is 
determined for each value of A a considered, and 
the optimal combination of shifts is identified by 
searching for the global minimum of the % 2 /dof. 
The ratios R a and Rz are related to A a and Az 
by R a = exp(-A a ); R z = i? 2 exp(-A z ). 

We considered matching the lattice data using 
q = 27m/ L as the momentum variable. The min- 
imum value for x 2 /dof of about 1.7 was obtained 
for R a ~ 0.815. This value for R a is considerably 
higher than the value of 0.716 ± 0.040 obtained 
from an analysis of the static quark potential in 
From this discrepancy, as well as the rela- 
tively high value for x 2 /dof, we may conclude 
that the gluon propagator, taken as a function 
of q, does not exhibit scaling behaviour for the 
values of (3 considered here. 



2 




Figure 1. \ 2 P er degree of freedom as a function 
of the ratio of lattice spacings for matching the 
small and fine lattice data, using q as the momen- 
tum variable. The dashed line indicates the ratio 
Rz of the renormalisation constants. 



Fig. [j] shows the result of the matching us- 
ing q = 2sin(g/2) as the momentum variable. 
This gives much more satisfactory values both 
for % 2 /dof and for R a . The minimum value for 
X 2 /dof of 0.6 is obtained for R a — 0.745. Taking a 
confidence interval where x 2 /dof < Xmin + ^ gives 
us an estimate of R a — 0.745137, which is fully 
compatible with the value [| of 0.716 ± 0.040. 



2. Model fits 

We have demonstrated scaling in our lattice 
data over the entire range of q 2 considered, and 
will now proceed with model fits. The following 
functional forms have been considered: 
Gribov §] 

L(q 2 ,M 1R ) 



W) = - 

Stingl § 

D r (q 2 ) 



q 4 + 2A 2 q 2 
Marenzoni [||] 

f 



L(q 2 ,M m ) 



D r {q z ) 



(q 2 ) 



2\l+a 



Ml 



Cornwall 



D r {q 2 ) 



(q 2 + M 2 ( q 2 ))ln 



q 2 + 4M 2 (q 2 ) 
"A 2 



(5) 



(6) 



(7) 



(8) 



Model A 

D r (q 2 ) = 
Model B 

D r (q 2 ) = 



A 



(<Z 2 +M2 R )i+« 
A 



D vv (q 2 



(q 2 y+° + (Miy+° 

Model C 

D r (q 2 )=Ae-^l M ^ +D vv (q 2 ) 
where D(q 2 ) = ZD r (q 2 ) and 



D vv (q 2 



D vv {q 2 ) = 



1 



q 2 + M 2 , 



-L(q 2 ,M uv ) 



1 



L(q 2 ,M)= ( -H(q 2 +M 2 )/M 2 ) 



-dr 



M(q 2 ) = M { In 



A 2 



4M 2 , , 4M 2 



-6/11 



(9) 
(10) 

(11) 

(12) 
(13) 
(14) 



We have also considered the models A and B, 
which are constructed from models A and B by 
setting M uv = M IR . Gribov's and Stingl's models 
(||) and @ are modified in order to exhibit the 
asymptotic behaviour of (g). Models A and B 
are constructed as generalisations of (0) with the 
correct dimension and asymptotic behaviour. 

All models are fitted to the large lattice data 
using the cylindrical cut defined in |l],f7j . The low- 
est momentum value was excluded, as the volume 
dependence of this point could not be assessed. In 
order to balance the sensitivity of the fit between 
the high- and low-momentum region, nearby data 
points within A(qa) < 0.05 were averaged. 

The x 2 P er degree of freedom and parameter 
values for fits to all these models are shown in 
table [l]. It is clear that model B accounts for the 
data better than any of the other models. The 
best fit to this model is illustrated in fig. |[ 

3. Discussion and outlook 

We have demonstrated that the gluon propaga- 
tor exhibits scaling over a wide range of momenta 
q. The data are consistent with a functional form 
D r (q 2 ) - D v 

1 



£> uv , where 



f 



2 q 4 + M 4 ' 



(15) 
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Table 1 



Parameter values for fits to models (||)-([lT]). The values quoted are for fits to the entire set of data. The 
errors denote the uncertainties in the last digit (s) of the parameter values which result from varying the 
fitting range ||. 



Model 


X Vdof 


Z 


A 




M uv or A 


a 


Gribov 


1972 


2.19±?| 




0.23l\ 4 






Stingl 


1998 


2.2 





0.23 






Marenzoni 


163 


2.41+° 2 




0.14+ 4 4 




O.29I2 


Cornwall 


50.3 


6.5±S 






0.271? 




A 


2.96 




1.2411, 


0.46+14 


0.961 4 ? 


i-silil 


A' 


3.73 


1.713 


0.84+2° 9 


0.48+! 2 7 




1.52li? 


B 


1.57 


1.78+g 


0.49te 7 


0.43l'i 


0.20l?9 


O.95I5 


B' 


4.00 


I.62I4 


0.581? 


O.40I2 




0.921Y 
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47 


2.0915° 


29+ 166 


0.22l 1 ° 6 


0.14+8 


0.49li° 6 
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Figure 2. The gluon propagator multiplied by q 2 , 
with nearby points averaged. The line illustrates 
our best fit to the form defined in ( |l0| ) . The scale 
is taken from the string tension S. 



M = 0.81_ 2 GeV, and D uv is the asymptotic 
form given by (|l2|) and (|l|). A more detailed 
analysis || of the asymptotic behaviour reveals 
that the one-loop formula @) remains insufficient 
to describe data in the region q 2 ~ 25GeV 2 . 

Issues for future study include the effect of Gri- 
bov copies and of dynamical fermions. We also 
hope to use improved actions to perform realis- 
tic simulations at larger lattice spacings. This 
would enable us to evaluate the gluon propagator 
on larger physical volumes, giving access to lower 
momentum values. 
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